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We study the expansion of a rotating, superfluid Fermi gas. The presence and absence of vortices 
in the rotating gas is used to distinguish superfluid and normal parts of the expanding cloud. We 
find that the superfluid pairs survive during the expansion until the density decreases below a 
critical value. Our observation of superfluid flow at this point extends the range where fermionic 
superfluidity has been studied to densities of 1.2 x 10 11 cm -3 , about an order of magnitude lower 
than any previous study. 
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Ultracold atomic gases have been used to create novel 
quantum many-body systems ranging from Bose-Einstein 
condensates, Mott insulators in optical lattices, to high- 
temperature supernuids of strongly interacting fermions. 
These systems offer a high degree of control over physical 
parameters including interaction strength and density. 
Many important features in these gases have a spatial 
scale too small to be resolved while the gas is trapped. 
A standard technique to reveal this physics is to switch 
off the confining potential and release the gas from the 
trap. A non-interacting gas expands ballistically and the 
expansion reveals its momentum distribution. The ex- 
pansion dynamics of an interacting gas is modified by the 
effect of collisions. This can result in classical hydrody- 
namic flow and in this case the expansion serves as a (not 
necessarily linear) magnifying glass for the trapped state. 
In contrast to classical hydrodynamics, superfluid hydro- 
dynamic flow does not rely on collisions. When a weakly 
interacting Bose-Einstein-Condensate (BEC) is released 
from an anisotropic trapping potential, superfluid hydro- 
dynamics leads to an inversion of the aspect ratio, often 
regarded as a hallmark of Bose-Einstein condensation Q . 

The expansion dynamics of strongly interacting Fermi 
gases has been the subject of a longstanding debate. For 
a weakly interacting ultracold Fermi gas anisotropic ex- 
pansion has been proposed as a probe for superfluidity, 
analogous to the case of weakly interacting BECs |2|. 
Anisotropic expansion has been experimentally observed 
in strongly interacting Fermi gases I n this case, 

however, the inversion of the aspect ratio can occur due 
to collisions between the expanding atoms even if they 
were initially at zero temperature jy, Q • So far exper- 
iments have not been able to discriminate between su- 
perfluid and collisional hydrodynamics, and indeed one 
would expect both effects to contribute: In the BCS- 
regime, the superfluid transition temperature Tc de- 
pends exponentially on the density. Starting at T < Tc, 
the superfluid gas should first expand according to super- 
fluid hydrodynamics. As the density drops, T approaches 
Tc and superfluidity cannot be maintained. From this 
point on, the gas should expand according to collisional 



hydrodynamics or enter a regime intermediate between 
collisional hydrodynamic and collisionless expansion. 

In this paper we study the expansion of a superfluid 
Fermi gas, in the regime where pairing is purely a many- 
body effect. We have observed superfluid flow even after 
5 ms of expansion, when the cloud size had increased by 
more than a factor of 4 and the peak density had dropped 
by a factor of 17 compared to the in- trap values. 

Superfluidity in Fermi gases has previously been estab- 
lished through the observation of vortex lattices HQ. To 
detect vortices in a rotating fermion pair condensate the 
pairs are transferred into stable molecules by sweeping 
an external magnetic field across a Feshbach resonance 
shortly after the gas is released from the trap 10]. Vor- 
tices can be observed only when the gas is still a super- 
fluid at the moment of the magnetic field sweep. At the 
final magnetic field (on the BEC side of the Feshbach 
resonance) the interactions are much weaker. Therefore 
the vortex core has higher contrast and is larger than 
near resonance. If the gas is no longer superfluid at the 
time of the field ramp, we expect the vortex core to fill in 
quickly and disappear. The observed vortex cores there- 
fore serve as markers for the regions which are superfluid 
at the time of the magnetic field ramp. 

Our experimental setup has been described earlier [l]| . 
Quantum degenerate fermionic 6 Li was prepared in an 
optical dipole trap after laser cooling and sympathetic 
cooling by 23 Na . An equal spin mixture of the two low- 
est 6 Li hyperfine states was created by a radio-frequency 
sweep at a magnetic bias field of 885 G. These states, 
labeled |1) and |2), exhibit a broad Feshbach resonance 
centered at a magnetic field Bo w 834 G. At magnetic 
fields below (above) £? , on the BEC (BCS) side, the 
scattering length a is positive (negative) and a nearby 
molecular bound state exists (does not exist). A fermion 
pair condensate containing about 5 x 10 6 fermion pairs 
was obtained by evaporatively cooling the spin mixture 
while ramping the magnetic field to 812 G. Note that at 
this magnetic field, the bond length of the molecular state 
is larger than the interatomic spacing, and the fermion 
pairs are bound by many-body effects. The final radial 
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FIG. 1: Superfluid expansion of a strongly interacting rotating Fermi gas. Shown are absorption images for different expansion 
times on the BCS-side of the Feshbach resonance at 910 G (0.0, 1.0, 2.0, 3.0, 3.5, 4.0, and 4.5 ms) and 960 G (0.0, 0.5, 1.0, 
1.5, 2.0, 2.5, and 3 ms), before the magnetic field was ramped to the BEC-side for further expansion. The vortices served as 
markers for the superfluid parts of the cloud. Superfluidity survived the expansion for several milliseconds and was gradually 
lost from the low density edges of the cloud towards its center. Compared to 910 G (a = —7200 an), superfluidity decayed 
faster at 960 G (a = —5000 ao) due to the reduced interaction strength. The total expansion time remained constant pjj. The 
field of view of each image is 1.2 mm x 1.2 mm. 



and axial trapping frequencies were uo r = 2ir x 120 Hz 
and ou a = 2tt x 23 Hz, respectively. To observe vortices 
as a probe of superfluid flow, the gas was set in rotation: 
two blue-detuned laser beams were rotated symmetrically 
around the cloud for 1 s at an angular frequency of 2tt x 80 
Hz We allowed 500 ms of equilibration before the 
magnetic field was ramped (in 500 ms) to several probe 
fields on the BCS side of the resonance. Finally, we stud- 
ied the expansion of the rotating superfluid: The gas was 
released from the optical trap and expanded at the probe 
field for a variable "BCS-expansion" time tecs- To trans- 
fer the remaining fermion pairs into stable molecules the 
magnetic field was then lowered in 400 /is to 680 G [l2(. 
Here, the cloud was given several milliseconds of "BEC- 
expansion". For absorption imaging the magnetic field 
was raised to 730 G in 500 /is before the last 2 ms of 
time-of- flight. For most of the data the total time-of- 
flight was chosen to be 11 ms 13]. An absorption image 
of the gas was obtained separately at tecs to determine 
the peak density and the peak Fermi momenta kp before 
the magnetic field sweep. 

Fig.[I]shows absorption images taken as outlined above 
for seven different BCS-expansion times at both 910 G 
and 960 G. The presence of vortices proves that super- 
fluid fermion pairs survived in the expanding gas for sev- 
eral milliseconds. As the density of the gas dropped dur- 
ing the BCS-expansion the vortices were gradually lost 
from the low density edges of the cloud towards its center. 
After 4.5 ms time-of-flight at 910 G and 3 ms at 960 G all 
of the vortices had decayed. If we regard the number of 
vortices as an indicator of the superfluid fraction of the 
gas, we can draw the "phase diagram" of Fig. O Here the 
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FIG. 2: (color online) "Phase diagram" of an expanding, ro- 
tating Fermi gas: At a given magnetic field the number of 
vortices served as a measure for the size of the superfluid re- 
gion in the gas. The number of vortices is plotted versus 
1/kF and l/\a\. The contour plot was created from a total of 
53 data points. In this diagram lines of constant kFd corre- 
spond to hyperbolas. The vortices decayed when the density 
(increasing 1/&f) or the scattering length (increasing l/\a\) 
was reduced. For weaker interactions, at smaller scattering 
lengths \a\, vortices were lost already at higher densities. The 
four data points shown mark the breakdown of superfluidity 
and are the same as those given in Fig. El (filled circles). 
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FIG. 3: The peak interaction strength at the end-points of 
superfluid expansion. Starting at a peak kFd in the optical 
trap (red triangles) , the vortices survived up to a critical peak 
kFd of -0.8 +/- 0.1 (filled circles), almost independent of the 
magnetic field (scattering length). The arrows mark the evo- 
lution during expansion. The critical kFd was obtained for 
each magnetic field by taking the average of the peak Jzf of the 
last cloud that still contained vortices at a time tecs and the 
peak kF of the cloud that did not contain vortices anymore 
500 {is later. The error in kFd is about 10% and dominated 
by the systematic error in the atom number. 



FIG. 4: Disappearance of bimodality. Zero temperature 
Thomas-Fermi profiles were fitted to the density profiles ob- 
tained after BCS-expansion at 885 G and subsequent BEC- 
expansion at 680 G. The % 2 of the fit was monitored as a 
function of the BCS-expansion time £bcs- A high \ 2 indi- 
cates a bimodal density distribution. Vortices were still ob- 
served after 5 ms of expansion (indicated by the dashed line in 
the figure) while the bimodality had already disappeared (for 
X 2 values smaller than 0.01 bimodality cannot be discerned). 
Hence, the absence of bimodality does not imply the absence 
of superfluidity. 



number of vortices is shown as a function of the inverse 
scattering length 1 /a and the inverse peak Fermi momen- 
tum 1/kp. As 1/kp increases at a given magnetic field, 
corresponding to the decrease in density during time-of- 
flight, vortices are lost. The reduction in the number 
of vortices for decreasing \a\ reflects the decrease of the 
superfluid fraction for smaller attractive interactions at 
a given temperature. In addition, the increase in the 
normal fraction leads to higher damping of the remain- 
ing vortex number 8]. Most importantly, however, we 
see that vortices are lost earlier in time-of-flight as the 
interactions are reduced. 

At all magnetic fields, we find that the peak interac- 
tion strength at the point where all vortices were lost 
is about constant, kpa ~ —0.8 (see Fig. |2j). As shown 
in Fig. [I] the loss of vortices occurred gradually and the 
surviving vortices were located within a circle of decreas- 
ing radius. We assume that the critical value of hp a for 
which superfluidity was lost, was first reached at the edge 
of the cloud and subsequently further inward. However, 
we were not able to confirm this picture quantitatively 
without a model that describes how the shape of the 
cloud and the bimodality develop during and after the 
magnetic field sweep. 

It is remarkable that the observation of superfluidity 
and fermion pair condensation for trapped gases has also 
been limited to values of kp\a\ larger than 1 pi. HfiL Ht^. 
This suggests that the underlying reason for this limita- 
tion is the same for a trapped and an expanding gas. One 
obvious scenario for the decay of the vortex lattice dur- 
ing expansion is the breakdown of superfluidity when a 
critical interaction strength is reached. As the density 
decreases, Tc/Tp drops while T/Tp remains constant 



(since the phase space density n x T -3 / 2 is invariant dur- 
ing expansion). Therefore Tc eventually becomes smaller 
than T everywhere in the cloud and superfluidity is lost. 
This scenario implies that the superfluid state evolves 
adiabatically during expansion, which is plausible: Even 
when the critical kpa is reached, the pair binding en- 
ergy still changes at a slower rate, A/ A, than the rate 
at which the pairs can respond to this change, A/ft 
Here A = (2/e) 7 ^ 3 Ep exp(— 7r/2kp\a\) is the pairing gap 
in the BCS limit (valid for kp\a\ ^ 1) where the 
peak Fermi energy Ep = h 2 kp 2 /2m and kp are density 
dependent. When superfluidity is lost, superfluid hydro- 
dynamics is probably replaced by an expansion interme- 
diate between collisional hydrodynamic and collisionless. 
For weakly interacting BECs, the decay of vortex lattices 
at finite temperature was studied theoretically in [2p| . 
and remarkably similar structures are found. 

Another explanation for the loss of vortices is a possi- 
ble failure of the transfer of correlated fermion pairs into 
molecules since the size of the fermion pairs increases 
with decreasing density. When the fermion pair size be- 
comes larger than the interparticle spacing, molecules 
might be formed out of uncorrelated nearest neighbors 
rather than out of correlated pairs. The magnetic field 
sweep then destroys the coherent many-body wavefunc- 
tion. 

Vortices || and bimodal density distributions 0, 
[Tjj have been used as indicators for superfluidity and 
pair condensation, respectively. If a fermion pair con- 
densate is transferred to the BEC side before its inter- 
action energy has been converted into kinetic energy, it 
continues to expand with the drastically reduced mean- 



FIG. 5: Loss of vortex contrast on resonance at 834 G. Shown are absorption images after a fixed total time-of-flight, but for 
different expansion times on resonance (2, 2.5, 3, 3.5, 4, 5, and 6 ms) before the magnetic field was swept to the BEC-side for 
further expansion. A gradual loss of the vortex contrast from about 15% (after 2 ms of expansion on resonance) to 3% (after 
5 ms) was observed across the whole cloud. The field of view of each image is 1.2 mm xl.2 mm. 



field energy of a molecular BEC at 680 G. This results in 
a clear separation of condensate and thermal cloud after 
further BEC-expansion. If the transfer of fermion pairs 
into molecules is delayed after releasing the gas from the 
trap, the fermion pair condensate initially expands just 
like the normal part of the cloud. This eventually leads to 
a loss of bimodality in the density profiles after the trans- 
fer. We can now study how the two indicators, vortices 
and bimodality, are related in this experiment. For short 
BCS-expansion tecs our data showed bimodality as well 
as vortices. However, the bimodality was gradually lost 
and could not be discerned after a longer BCS-expansion 
although vortices were still visible (see Fig. El for details). 
The absence of bimodality therefore does not indicate a 
breakdown of superfluidity. 

So far we have studied the expansion of the gas on the 
BCS side of the Feshbach resonance. On the BEC side 
and on resonance, Tc is proportional to Tp so that T /Tc 
is constant during expansion. Therefore, one would not 
expect to observe a breakdown of superfluidity in expan- 
sion. Fig. shows absorption images that were obtained 
after an initial expansion of the cloud on resonance at 
834 G. In contrast to the situation on the BCS-side of 
the resonance no vortices were lost. Instead, the vortex 
contrast decreased uniformly across the cloud for longer 
expansion times. Vortices have been detected at total 
densities as low as 1.2 x 10 11 cm -3 in the wings of the 
expanded cloud. Here the critical temperature Tc of ap- 
proximately 0.2 T F was below 20 nK (k B T F is 
the local Fermi energy). We believe that the reduction 
in the vortex contrast is due to the low density of the gas 
after long BCS-expansion: after the magnetic field sweep 
the vortex cores cannot adjust quickly enough to the high 
contrast and large size they would have in equilibrium on 
the BEC-side. This loss of contrast limited our study of 
the breakdown of superfluidity to magnetic fields above 
880 G. 

In conclusion we have shown that superfluid pairs can 
survive during the expansion of a strongly interacting 
Fermi gas. This is the first observation of non-equilibrium 
superfluid flow in such systems. It has allowed us to 
observe fermionic superfluidity at total densities as low as 



1.2 x 10 11 cm -3 . An intriguing question for future studies 
is whether fermion pairs expanding from two clouds can 
coherently interfere. 
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